Nucleon participants or quark participants? 
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We show that centrahty dependence of charged particle pseudorapidity density at midrapidity 
in Au+Au collisions at RHIC is well described as proportional to the number of participating 
constituent quarks. In this approach there is no need for an additional contribution from hard 
processes usually considered in the models based on the number of the nucleon participants. 



o 
o 

(N 



> 

O 
(N 
O 
cn 
o 



X 



PACS numbers: 25.75.-q, 25.75.Ld 

Keywords: centrality, participants, constituent quark 



I. INTRODUCTION 

Charged particle multiplicity densities near mid rapid- 
ity in high energy nuclear collisions depend strongly on 
collision centrality. In order to better understand this 
dependence and to disentangle pure nuclear effects, this 
density is often plotted as per nucleon participant pair . 
Participants are the nucleons that have encountered at 
least one inelastic collision. The number of participants 
at a given centrality is usually calculated in the Glauber 
model either analytically or with the help of a Monte- 
Carlo algorithm. The charged particle density per par- 
ticipant increases with centrality. In the 5% most cen- 
tral Au-f Au collisions it is about 20% larger than that 
in semi-peripheral collisions (50-70% centrality region), 
and it is about 50% larger compared to pp collisions at 
the same energy. 

The reason for this increase in number of produced 
particles per participant at midrapidity is still not well 
understood. The most common explanation of the phe- 
nomena involves particle production in hard processes. 
Hard processes have much smaller cross-sections than 
soft collisions and depend differently on collision central- 
ity. They scale with the number of binary collisions (the 
number of collisions the incident nucleon would experi- 
ence if it were not altered at all while passing through the 
nucleus). The number of binary collisions increases with 
centrality faster than the number of participants; this re- 
sults in an increase of particle production per participant 
nucleon as centrality increases. In such approaches the 
particle density is often presented simply as 01: 



dN^h 
dy 



(X aN^-part + (1 - a)NN-coll , 



(1) 



where the parameter a is the relative fraction of particles 
produced in the soft collisions, and (1 — a) is the relative 
fraction produced in hard collisions. With proper pa- 
rameters, this fits the data fairly well; see Ref. 0]. Note, 
however, that in such models, the relative contribution of 
hard processes is expected to increase with the collision 
energy. The data seem inconsistent with such an energy 
dependence. 

In the approach proposed in this paper, both nuclei and 
single nucleons are considered as a superposition of con- 



stituent quarks (also often called as "dressed" quarks or 
valons); there are three such dressed quarks per nucleon. 
The concept of the constituent quarks has been known for 
many years (see Refs. and references therein). The 

constituent quark approach is able to explain many fea- 
tures of hadron-hadron and hadron- nucleus collisions |^ . 
QCD calculations support the statement that inside a nu- 
cleon there are three objects of the size of 0.1 - 0.3 fm (see 
Ref. and references therein) ; for some recent works us- 
ing or discussing constituent quarks, see also Refs. 

In the constituent quark picture, a NN collision looks 
like a collision of two light nuclei. Most often only one qq 
pair interacts, with other quarks being spectators. Only 
part of the entire nucleon energy is spent for particle pro- 
duction at midrapidity (as y/Sqq ~ v^saw/3). The quark 
spectators form hadrons in the nucleon fragmentation re- 
gion. In the case of AA collisions, more than one quark 
per nucleon interacts due to the large nucleus size and 
the possibility for quarks from the same projectile nu- 
cleon to interact with different target nucleons. The goal 
of this study is to find the number of produced particles 
per participant quark (pair) and to check for its centrality 
dependence. 



II. CALCULATIONS OF THE NUMBER OF 
PARTICIPANTS 

We calculate mean number of nucleon/quark partici- 
pants using a Monte-Carlo based implementation 9] of 
the nuclear overlap model 0. We use Woods-Saxon 
nuclear density profile 



UAir) = 



no 



l + exp[(r-i?)/d]' 



(2) 



0.17 fm~^i? 



(1.12^1/3 



with parameters: uq — 
0.86^/3) fm, d = 0.54 fm. 

In the nuclear overlap model, the mean number of par- 
ticipants in j4 + _B collision at impact parameter b is given 
by 



part,. 



AB = / d^S Ta{3){1 - [1 



'b ^ ^ 
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dhTB{s){l-[l- 



aNNTAjs- b) 
A ' 



(3) 



where T{b) ~ J°°^ dz UAi^JlP' + z^) is the thickness func- 
tion; then (1 — (TjvAfT'yi(^)/^)'^ is the probabihty for a nu- 
cleon to pass through the nucleus without any coUision. 
We use the inelastic cross section u^n = 41 mb at 
./sWn = 130 GeV. 

The number of participating nucleons for a given cen- 
trality can be determined directly using the web inter- 
face . In order to calculate the number of participating 
quarks we downloaded FORTRAN code and modified 
it by increasing the density three times [n\ — 3no = 
51 fm~'^) and changing a^N to Oqq. For the quark- 
quark cross-section in our calculation, we use two values 
Oqq — unn 1^ — 4.56 mb and a somewhat arbitrary value 
of (Jqq — 6 mb; the latter was used mostly to illustrate 
the sensitivity of the results to the value of aqq. The 
choice of aqq = 6 mb is not unreasonable since at RHIC 
energies, approximately 1.2 - 1.3 quarks per nucleon can 
participate in a single A^A^-coUision 5]. In principle, aqq 
could be probably as high as 8 mb based on the early 
estimates of {vq/RN)^ ~ 1/5 4]. 

Fig. ^ shows the number of the nucleon and quark 
participants vs. impact parameter. Fig. ^ presents the 
centrality dependence of the ratio of A'g-part/A^Af-part. 
Smooth curves are the polynomial fits to the Monte-Carlo 
results to smooth out the statistical fluctuations. The ra- 
tio Nq-part/N^-part IS used later for the renormalization 
of the particle pseudorapidity distributions from per nu- 
cleon participant to per quark participant. 
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FIG. 1: Impact parameter dependence of (a) the number of 
the nucleon and the quark participants, and (b) the ratio of 
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Fig. 121 presents the comparison of our calculation of 
^N-part in the nuclear overlap model with PHOBOS cal- 
culations based on the HIJING model. The number of 
participants is plotted as a function of a given fraction of 
the total cross section. Open circles represent PHOBOS 
calculations. The nuclear overlap model results (using 
the Woods-Saxon density profile, the same as used in 
the HIJING model) are shown by solid symbols. Note 
a small deviation of our calculations from that of PHO- 
BOS in very central region; this is not important for the 
conclusion of the current study. 
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FIG. 2: Mean number of nucleon participants vs centrality in 
the nuclear overlap model (solid symbols) and from PHOBOS 
calculations (open circles). 

The PHOBOS Collaboration presents their results on 
centrality dependence of the charged particle pseudora- 
pidity density by plotting it vs. the number of the nu- 
cleon participant pairs. In this paper, we continue to use 
the same quantity for the centrality characterization, but 
note that in the constituent quark picture, the number 
of the nucleon participants no longer has the meaning of 
the number of the particle production sources. 



III. RESULTS 

The centrality dependence of the charged particle mul- 
tiplicity per participant pair is shown in Fig. |3| The re- 
sults per nucleon participant pair are in the upper part 
of the figure and the results per quark participant pair 
are shown in the lower part. The original PHOBOS re- 
sults on dNch/drj per nucleon participant pair (calcu- 
lated using HIJING) for = 130 GeV and 200 GeV 
are shown in solid symbols. In open symbols we also 
show the same results renormalized for the number of 
the nucleon participant pairs from our calculations using 
the nuclear overlap model. Our main result - the cen- 
trality dependence of the charged particle pseudorapidity 
density per quark participant pair ~ is presented in the 
same plot. 
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FIG. 3: Nch per nucleon and quark participant pair vs cen- 
trality. The results for quark participant pair are shown for 
aqq = 4.56 mb (solid symbols) and aqq = 6 mb (open sym- 
bols). 



We observe no, or even slightly decreasing, dependence 
of {dNch/dy)/Nq-part on centrality, with the ratio being 
dependent only on the energy of the collision. The slight 
decrease in particle production at midrapidity with in- 
creasing centrality could be due either to low values of 
the constituent quark inelastic cross section used in our 
calculations or to parton saturation effects. Note that 
in the constituent quark picture, {dNch/dy) /Nq-part as 
a function of centrality depends very weakly on the col- 
lision energy, as the change in the inelastic cross section 
is probably less than 5% between ^Jsnn = 130 GeV and 
200 GeV. 

Hard processes scale with the number of binary col- 
lisions. Although it was not necessary to include the 
contribution of hard processes into our calculation in or- 
der to describe the centrality dependence of the charged 
particle density at midrapidity, we have calculated the 
number of binary collisions as well: see Fig. 0] Note 



that the number of binary collisions per participant has 
a much weaker centrality dependence in the constituent 
quark approach than it has in the the nucleon participant 
model. 



IV. SUMMARY 

We have shown that the particle multiplicity density 
near mid rapidity scales linearly with the number of con- 
stituent quark participant pairs. The experimentally ob- 
served increase of dNch/drj per nucleon participant pair 
with centrality in this picture is explained by the relative 
increase in the number of interacting constituent quarks 
in more central collisions. 
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FIG. 4: Centrality dependence of the ratio of binary collisions 
to the number of nucleon and quark (upper curve corresponds 
to aqq — 6 mb and lower to aqq = 4.56 mb) participants. 
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